Experimental imaging and atomistic modeling of electron and hole quasiparticle wave functions in 

InAs/GaAs quantum dots. 
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We present experimental magnetotunneling results and atomistic pseudopotential calculations of quasiparticle 
electron and hole wave functions of self-assembled InAs/GaAs quantum dots. The combination of a predictive 
theory along with the experimental results allows us to gain direct insight into the quantum states. We monitor 
the effects of (i) correlations, (ii) atomistic symmetry and (iii) piezoelectricity on the confined carriers and (iv) 
observe a peculiar charging sequence of holes that violates the Aufbau principle. 



I. INTRODUCTION 

The localized, quantized and entangled states of carriers 
confined to quantum dots manifest a wealth of novel physical 
phenomena that are observed experimentally mostly through 
measurements of the energies of characteristic processes such 
as the formation, splitting, and charging of excitonsif^iM. 
However, yet another level of control might be achieved 
through engineering of the wave functions themselves through 
manipulation of their degree of localization, spatial anisotropy 
or angular-momentum character The first crucial capability 
that would allow the design of a device on the fundamental 
level of its wave function character is wave function moni- 
toring along with an understanding of the controlling physical 
parameters. However, both experimental imagin g^'^'^'^'^'^" and 
theoretical calculations of the many-particle wave functions 
of dots present a formidable challenge. From the theoretical 
point of view, the tradition has largely been to fit a few mea- 
sured energy levels by adjusting a few parameters in the con- 
fining potential within simple models (e.g parabolic effective- 
mass models), or directly adjusting the energy-related quan- 
tities (e.g., tunneling in dot molecules or fine-structure split- 
tings), without imposing physical reality on the wave func- 
tions other than their boundary conditions. However, theoreti- 
cal determination of wave functions are more challenging than 
modeling of the corresponding eigenenergies because of their 
high sensitivity to subtle physical effects. For instance, piezo- 
electric terms or particle-particle correlations have a rather 
small effect on energies" yet can affect wave function shapes 
rather clearly. In this contribution, we present a combina- 
tion of experimental magnetotunneling technique and a pre- 
dictive theoretical modeling of wave functions that gives us 
an unprecedented insight into the physics of carriers in con- 
fined geometries. This combined approach provides a direct 
answer to important physical problems such as the hole fill- 
ing sequence'^ '^ that violates the Aufbau principle, the im- 
portance of the atomistic symmetry, correlations and piezo- 
electricity and we address the decade-old question'^ whether 
the anisotropy observed in the spectroscopy of self-assembled 
dots'^ is caused by shape anisotropy (oval dots) or piezoelec- 
tric fields. We believe that the type of studies descried here 



are very promising for wave function engineering and design 
of future nano- scale devices. 



II. EXPERIMENTAL METHOD 

The systems studied here are self-assembled InAs quan- 
tum dots, grown by the Stranskii-Krastanov method'''-. 
They are investigated by capacitance-voltage (C-V) 
spectroscopy'^'-'—. The energy of the quantum dot 
states are shifted with respect to a carrier reservoir (back 
contact) by the applied voltage. Characteristic maxima in the 
capacitance appear, each time an additional electron (or hole) 
can tunnel into the dots. This makes it possible to determine 
the addition spectrum of such dots with great accuracy'^. 

To map out the probability densities of the quantum dot 
states, we extend the approach of Vdovin et al^ and mea- 
sure the ac-tunneling current between the dots and and the 
back contact. Additionally, a magnetic field B is applied per- 
pendicularly to the tunneling direction. This field imposes an 
in-plane momentum 
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on the tunneling of the carriers, where d is the tunneling dis- 
tance. At sufficiently high frequencies, the amplitude of a 
capacitance maximum is a measure of the tunneling rate be- 
tween the dots and the back contact'^. The tunneling rate, in 
turn, is proportional to the probability density of the quasi- 
particle wave function in k-space^''''^". Experimentally, map- 
ping out the wave functions thus requires recording the C-V 
amplitudes as a function of the in-plane magnetic field for dif- 
ferent azimuthal orientations^. 

The investigated samples are GaAs-(Alj:Gai_j;As) Schot- 
tky diodes, grown by molecular beam epitaxy, with embed- 
ded InAs quantum dots. The p-doped (hole) sample was pre- 
pared as described in Ref. [13 however, with a slightly thicker 
tunneling barrier of 19 nm, to facilitate wave-function map- 
ping. The layer sequence of the n-doped sample for the elec- 
tron spectroscopy can be found in Wibbelhoff et at.—. Schot- 
tky diodes were prepared by alloying ohmic contacts and de- 
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positing Cr-Au topgates (300/im x 300 fim). The C-V spec- 
troscopy was carried out using a standard LCR meter (Agilent 
4284A) with an ac voltage modulation of AV — 10 mV. The 
frequency was appropriately chosen (8 - 40 kHz) so that the 
capacitance amplitude reflects the tunneling rate^ '^. To de- 
termine the tunneling probability as a function of fc| | , and this 
way map out the quasi-particle wave function in momentum 
space, C-V spectra for in-plane fields up to B = 26 T and 
for azimuthal angles in steps of 15° (starting parallel to [Oil] 
crystal direction) were evaluated. The in-plane momentum 
follows from Eq. ([T]i. The normalized C-V amplitudes of the 
different charging peaks (0 to 6 carriers per dot) are plotted as 
a function of fc|| in Figs. [Tta) and (c). 



4'/5(xi, ■ • ■ , xat) is a Slater determinant of N electrons, and 
Cp^'^ is its weight. Therefore, 

M,,N ^Y.'^N ~ l\c,\N){^y,\^,) 

i 

= EE^^''"^^r(*i"'"^i'5^i<^) (5) 

i a, (3 



with 



1 if|ci>(^-i))=c,:|$W), 



(*r-^)|fid*r> = <l -1 if|a>(^-^') = -c.|<'), 



otherwise. 



III. THEORETICAL METHOD 

A full theory of magneto-tunneling would require a self- 
consistent calculation of the transport properties of the full 
system under external field. This is still prohibitive at an atom- 
istic level for such large systems (~ 10® atoms). Instead, we 
used a simplified transmission theory, which ignores the non- 
linear effects of electric field and device structure effects and 
assumes resonant tunneling, i.e., the emitter state are tuned in 
such a way to be in resonance with the quasiparticle quantum 
dot state. We have calculated the transition rate of an electron 
or hole from an emitter in state k to a quantum dot containing 
N particles following the work of Bardeen^'. In this approx- 
imation the transition rate is given by T^.a? |A^K,Afp and 
the transition matrix elements for the transfer of one particle 
from the emitter in state k to the quantum dot state | A^), filled 
by electrons, is given by: 



(2) 



Here, (/'*(r) is the probing or emitter wave function and is 
generated by the external source, ^'qd (x) is the quasi-particle 
excitation between the A^ — 1 particle states | A^ — 1) and the 
A^ particle states | A^), i.e.. 



*QD(x)-^(Ar-l|c,|iV) V'.(x) 



(3) 



ipi (x) is the atomistic ith single-particle wave function and Ci 
is an electron (hole) annihilation operator. 

To obtain the correlated many-body states, we use the con- 
figuration interaction approach where the many-body wave 
function is written as a superposition of different Slater de- 
terminants (configurations), such as, 

a 

1^) = E^r^*/3(^i'--- '^^^-I'^jv), (4) 

where $a(xi, ■ • • ,XAr_i) is a Slater determinant of A^ — 1 
electrons, and C^~^'^ is its weight. Accordingly, 



where i is the single particle states. 

In the past, AA^.n was calculated by using single-band ef- 
fective mass wave functions for ip^. However, an effective 
mass theory ignores the atomistic character of the wave func- 
tions and may miss inter-band and inter-valley effects. Here 
we use an atomistic empirical pseudopotential approach that 
takes multi-band, multi-valley and spin-orbit couphng into ac- 
count. 

Our atomistic wave functions for the quantum dot states can 
be written as: 



Nb Nk 

^■'W = EE^n,U".k(x)e'' 



(6) 



where u„,k(x) are the strained InAs bulk Bloch wave func- 
tions (linear combination of bulk bands methods^^), Nb and 
Nk are the number of bands and k-points, respectively. The 
probing emitter wave function can be written as: 



(7) 



where ?2k(x) is the Bloch part of the emitter wave function, 
and is not precisely known. The projection in Eq. (|5]) can be 
written as: 



Nb 

E(^k|Un,k)ci'^k 

n 



(8) 



Since we do not know the exact form of the Bloch part of the 
probing wave function, we assume that (uklwn.k) — const. 



IV. IMAGING WAVE FUNCTIONS ESTABLISHES 
CHARGING SEQUENCE 

The two issues here are: (i) Whereas the single-particle or- 
bital energies follow the order S, P, D, the addition of carriers 
may not successively fill the levels in that oder, but skip one 
shell, violating the Aufbau principle. Furthermore, (ii) the P 
states may split into PI and P2, even if the geometric shape 
of the dot is perfectly cylindrical. This "symmetry breaking" 
results from the fact that even in perfectly cylindrical/lens- 
shaped dots made of zincblende material the atomistic sym- 
metry and the strain symmetry (driving piezoelectricity) is 
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C2v where PI and P2 need not be degenerate. Figures |2^) 
and b) show two different filling sequences for electrons that 
may result from these issues and that can be distinguished by 
wave function imaging. Indeed, the fourth electron in Fig.|2fa) 
tunnels into the PI state while it tunnels into the P2 state for 
filling "Sequence 11". Since PI and P2 have a different ge- 
ometrical shapes both scenario can be discriminated by map- 
ping out the tunneling amplitudes. 



V. EFFECT OF THE QUANTUM DOT SHAPE 

One of the challenges posed at the onset of any comparison 
between theory and experiment is given by the experimen- 
tal determination of the quantum dot shape, that will subse- 
quently be used in the simulations. According to AFM- and 
SEM-measurement of the uncapped dots, the shape seems cir- 
cular. However, it is well known that the shape of the dots 
changes significantly by the overgrowth process and during 
this process an anisotropy might be introduced as well as a 
smoothing of the interfaces by diffusion and exchange be- 
tween In and Ga. Note that the measurements are performed 
on a statistical ensemble of dots of slightly different shapes 
and sizes. 

For the calculations to be representative and to assess the 
robustness of our theoretical results, we surveyed a large num- 
ber of dot shapes (4 different heights, 3 different base sizes 
and 5 different elongations) and the qualitative features of the 
results remained. In the present contribution we will show re- 
sults for three different dots detailed in Table |T] and labeled 
as Dl, D2 and D3. Dl and D3 both have a circular base but 
different heights and D2 has an elliptical base with elongation 
along the [110] direction. From the magnetotunneling results, 
we have a hint for an elongation of the dot: In Figure |3] we 
show contour plots of the transition probability measured by 
magnetotunneling spectroscopy [(a) and (d)] and calculated 
for dot Dl [(b) and (e)] and dot D2 [(c) and (f)]. We show 
results for the tunneling of the first electron [(a),(b),(c)] and 
the first hole [(d),(e),(f)] into quantum dot states with strongly 
dominant orbital S'-character The experimental results show 
an elongation of the signal along the [110] direction in real 
space. The theoretical results for the elongated dot D2 agree 
well with the experiment while the dot Dl with a circular 
shape differs from the experimental picture. We interpret this 
result as a strong hint for a structural elongation and conclude 
that dot D2 is the one closest to the experimental reality. 



VI. ADDITION ENERGIES 

In Table [n] we summarize the experimental and calculated 
addition energies A(A^ — 1, N). We calculated the addition 
energies from our many-body energies resulting from config- 
uration interaction. We define A(A^ — 1, N) as the difference 
between the charging energies for N and {N — 1) particles. 



/i(iV) and /i(iV — 1) respectively: 

A{N-1,N) = fi{N) - fi{N - 1) (9) 
= E{N) - 2E{N - 1) + E{N - 2) (10) 

The charging energy fJ.{N) is the energy required to add an 
additional carrier to a dot already occupied by N charges. 

Both electron and hole addition energies in Table |ll] agree 
fairly well with the experimental results. As a general trend, 
the addition energies for the holes are underestimated by the 
theory while no such general underestimation can be observed 
for the electron addition energies. The difference in the bare 
magnitude of the addition energies can therefore not be at- 
tributed to a simple difference in size between the experimen- 
tal and the simulated dots. 

The failure of the addition energies to unambiguously iden- 
tify the best choice between dot Dl, D2 and D3 (table 
highlights how complementary the magnetotunneling results 
are in pinpointing salient features of the quantum dots, such 
as elongation not revealed in charging or addition energies. 



VII. RESULTS FOR ELECTRONS 

Figure [TJa) shows the experimentally determined quasi- 
particle wave-functions of electrons for different occupation 
numbers N. The wave functions of the lowest two states (N 
= 1 and N = 2) are clearly S'-like and their similarity suggests 
weak correlations in the N = 2 case. A slight elongation of 
the overall shape of the = 1, 2 states is discernible along 
the [110] direction in fc-space, corresponding to an elongation 
along [110] in real space, as discussed before. The next set of 
higher lying states = 3, 4 exhibit a node along [110], while 
the states for = 5, 6 are oriented perpendicularly, along 
[110]. Note that the experimental resolution only allows us 
to map the 3, 4-states and the 5, 6-states together, hence the 
states 3, 4 and the states 5, 6 are simply doubled in Fig.fTla). 

Figure [Ttb) shows our theoretical results for electrons in 
dot Dl (not elongated). We calculated the electron states 
for the elongated dot D2 as well and the results are qualita- 
tively the same (not shown). The result of Fig. [TJb) are in 
good agreement with the experimental findings. The theoret- 
ical and experimental results conclusively point to the filling 
"Sequence I" of Fig. |3a). This implies that the P1-P2 split- 
ting is large enough so that Pi is fully occupied before P2, 
following Aufbau. Had the P1-P2 splitting been smaller than 
the sum of exchange energy and the difference in Coulomb 
energies for one electron in Pi and another in P2, we would 
have obtained the charging sequence of Fig.|2b)^^ following 
the Hund rule that was observed in large dots^'*. We choose 
to show the results of the dot with circular base Dl (while 
dot D2 is closer to the experimental situation) to emphasize 
the fact that shape elongation is not necessary to explain the 
results and to warn about the tempting conclusion that split- 
ting of P-levels or wave function anisotropy is an indication 
of dot asymmetry. Our work reveals the importance of atom- 
istic symmetry (not shape): Since our theoretical results from 
Fig- Hfb) for dots with cylindrical base agree very well with 
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our experimental findings, there is no need to assume shape 
asymmetry. Thus, the orientation of the P-states in Fig.fTfb) 
is a resuh of the atomistic nature of the underlying zincblende 
crystal lattice in contradiction with effective mass models that 
lead to degenerate and isotropic P-states. 



Vm. RESULTS FOR HOLES 

Figure[TJc) shows the experimentally determined probabil- 
ities for the first six hole states. The data for = 1, 2 shows 
again the shape of an S-state with, as was the case for the elec- 
trons, a slight geometric elongation along [110] in reciprocal 
space suggesting a slight elongation of the dot along the [110] 
direction. Close inspection of the = 3 and the = 4 hole 
states shows that they both exhibit a node along the [110]- 
direction, whereas the = 5 hole state is almost circularly 
symmetric with a clearly developed minimum in the center 
and hole state A^ = 6 exhibits nodes along [110]. 

For the theoretical calculations of Fig. [Ifd), we used the 
elongated dot D2 (see TablelJi. 

The transitions of the first two holes Oh Ih, and Ih 
2h have no nodes and resemble the case of electrons. The gen- 
eral feature that the theoretical results show narrower, sharper 
peaks can be traced back to the fact that the experiment probes 
an ensemble of many quantum dots with slightly different 
shapes and hence slightly different transition energies, while 
the theory is performed assuming a single quantum dot. 

For the third and fourth holes the 2h 3h, and 3h 4h 
transition amplitudes are anisotropic with peaks developing 
along the [110] direction. This is the signature of the first hole 
P state. Our work reveals the importance of inter-particle 
correlation effects: Indeed, a closer analysis of the theoret- 
ical results for the 2h^3h transition shows that 91% of the 
initial state is given by the /iq -configuration (configuration as 
described in eq. (|4|i, i.e, both holes occupy the first single- 
particle hole level ho. 88% of the final state is given by the 
/iQft,}-configuration where two holes are in state ho with S 
orbital character and one hole is in hi with P orbital charac- 
ter. This analysis shows that the tunneling hole is of orbital 
P-character It is interesting to note how the mainly-single 
band electron and the multi-band hole P-states differ in their 
quasiparticle tunneling amplitude. 

For the fifth hole, 4h 5h is different and shows mostly 
isotropic features. This is the signature of the tunneling 
into the D state. Indeed, 85% of the initial state is given 
by the -configuration and 82% of the final state by the 
/ipftf/ig -configuration (where the last hole is in the D state 
/13). This shows, that the tunneling hole is mainly of orbital 
D-character Comparison of the data in Fig. [Tfc) and (d) 
shows that even subtle differences in the shape of the wave 
functions can be resolved experimentally: The qualitative dif- 
ference in the calculation between the 4h 5h and the 5h 
6h transitions (namely equal or different amplitudes along the 
[110] and [110] directions) are clearly reflected in the spec- 
troscopic data. To emphasize the significance and clarity of 
the signature we obtain experimentally and theoretically, we 
artificially simulated the tunneling into a P2 state instead of 



a D state in Fig. |4] The left panel shows the artificial sit- 
uation where we fixed the final configuration to a hyi\h\- 
configuration, where the last hole is in the P2 state /12, fol- 
lowing the Aufbau principle. The right panel is the repetition 
of our actual result where the hole tunnels into a D state. Both 
figures being very different, we conclude that the signature is 
a strong indication that, indeed, the _D-state is filled before the 
P2 state. This is one of our main findings. 

Recent theoretical calculations^^, showed, that even for a 
conventional filling of the shells, following the Aufbau prin- 
ciple, the magnetic field dependent charging may describe the 
experimental results of Reuter et al}^ therefore challenging 
the interpretation of an unusual shell filling. We believe that 
our magnetotunneling spectroscopy results give additional ev- 
idence to supports the scenario of Reuter et ali^ and He et 
al. of a violation of the Aufbau principle. 

For the sixth hole, 5h 6h, regains some anisotropic char- 
acter with stronger maxima along the [1 10]-directions. The fi- 
nal state is still mainly given by holes in D states: h^h\h\, but 
now to a somewhat lower percentage of 77%. The remaining 
23% are configurations that include P2 state that have max- 
ima along the [1 10]-direction. This is hence an effect of corre- 
lation that tend to become more important for heavily charged 
states. Our work reveals the importance of piezoelectricity 
on the shapes of wave functions: While it only weakly af- 
fects the energies" it modifies the single-particle wave func- 
tions significantly. Ignoring the piezoelectric effect leads to 
an anisotropic 4h^5h transition and a more isotropic 5h^6h 
transition, in contrast to the correct theoretical treatment and 
the experiment. 



IX. ROBUSTNESS OF THE RESULTS 

To illustrate the robustness of our results we depict in Fig.|5] 
the theoretical results we obtain for hole charging on dot Di 
that possesses a circular base. We believe that the slightly 
elongated dot D2 from Fig. [Ttd) fits better the experimental 
situation, as discussed in section [V] but we want to illustrate 
the effect of shape on the results. The results in Fig.|5]are qual- 
itatively very similai" to Fig. [TJd) for Oh lh,lh 2h,2h 

3h and 3h 4h. For the ti'ansition 4h 5h the results 
are less isotropic in DZ than for the elongated dot D2 and 
shows more pronounced maxima along the [110] direction. 
The comparison of this transition with a tunneling of the hole 
into the second P-state (left panel Fig. |4|i show an even more 
pronounced difference than for the elongated dot D2: the ori- 
entation of the peaks are rotated by 90°. Hence, the signature 
of tunneling into an orbital Z?-state is even stronger in a dot 
with circular base than for the elongated dot D2. However, 
the experimental results agree better with the calculations of 
dot D2 than 1)3, as we already concluded in section [V] The 
results for 5h 6h are very similar for dot D2 and Di. 
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X. SUMMARY 

In conclusion, we measured and calculated the quasiparti- 
cle transition amplitudes resulting from a magnetotunneling 
experiment. Experiment and theory represent a novelty with 
the first experimental mapping of holes and the first atomistic 
calculation of quasiparticle transition amplitudes in quantum 
dots. We provide a direct evidence for the violation of the 
Aufbau principle for holes. The excellent agreement between 
the experiment and the theory allows us to further analyze 
and fully understand the nature of each of the many-body 
states probed. This analysis reveals the orbital character of 



the states, and shows the importance of the atomistic symme- 
try, correlations and piezoelectricity in the results. We believe 
that this type of approach delivers unprecedented insight into 
the electronic structure of self-assembled quantum dots and 
might be used for wave function engineering. 
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U.S. Department of Energy, Office of Science, Basic Energy 
Sciences, under Contract No. DE-AC36-99GO10337 (LAB- 
17) to NREL, from the DFG through the GRK384 and through 
the German Federal Ministry of Education and Research and 
under the "nanoQUIT" Program. 
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TABLE I: Label and dimension, in nm, for the 3 dots considered. All 
dots are pure InAs with an overall lens shape. For dot D2 the base is 
elliptical while dots Dl and Z)3 have a circular base. 



Dot 


base along 


[110] base along [110] height 


Dl 


25.0 


25.0 5.0 


D2 


26.0 


24.0 3.5 


D3 


25.0 


25.0 3.5 



TABLE IL Addition energies for electrons and holes in meV. The ex- 
perimental values of electron and hole addition energies are extracted 
from Refs. [Tsl and[T2l respectively. "Theory Dl" are the results for 
an InAs lens-shaped dot with circular base of 25 nm diameter and 
5 nm height. "Theory D2" is for an InAs lens-shaped dot with an 
ellipsoidal base of 26 nm (along[110]]) x 24 nm (along [110]) and 
3.5 nm height. "Theory D3" is for an InAs lens-shaped dot with a 
circular base of 25 nm diameter and 3.5 nm height. 



Addition 






Electrons 




Energy 


A(l,2) A(2,3) A(3,4) A(4,5) A(5,6) 


Experiment 


22 


57 


11 


21 12 


Theory Dl 


22 


63 


19 


22 19 


Theory D2 


20 


64 


17 


15 17 


Theory D3 


21 


66 


16 


18 16 




Holes 


Experiment 


24 


34 


17 


23 15 


Theory Dl 


18 


21 


16 


21 14 


Theory D2 


11 


15 


13 


15 13 


Theory D3 


19 


23 


16 


22 14 
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FIG. 1: Figure can be found at http;//www.sst.nrel.gov/nano_pub/mtS-prep rint.pd^ Color online) Experimental [(a) (c)] and theoretical 
[(b) (d)] quasi-particle probability densities for electrons [(a) (b)] and holes [(c) (d)]. The directions given by arrows and labeled as [110] and 
[1-10], are real-space crystallographic directions (while the plots are in k-space). The Calculations are performed on a single quantum dot 
while the experiment probes an ensemble of quantum dots. For the electrons (holes), the reciprocal lattice vectors span a range of -6 to -1-6 10* 
m~^ (-7 to +1 10* m"^). 



FTG. 2: Figure can be found at 

http://www.sst.nrel.gov/nano_pub/mts_preprint.pdf Single- 
particle orbital filling sequence. Two possible scenarios for electrons 
are given in [(a) (b)] and our results for holes in [(c)]. 



FIG. 3: Figure can be found at 

rhttp://www. sst.nrel.gov/nano_pub/mts_ preprint.pdf Color on- 
line) Contour plot of the transition probability for the transition Oh 
Ih. The directions are given in real space while the plots are in 
reciprocal space, i.e., dot D2 is elongated along the [110]-direction 
and so is the transition probability map. The reciprocal lattice 
vectors span a range of -5.4 to -1-5.4 10* m~^ for a) and of -4.3 to 
-1-4.3 10* m"^ forb) and c). 



FIG. 4: Figure can be found at 

pttp://www.sst.nrel.gov/nano_pub/mts_preprint.pdf 'Color On- 
line) (Left) Quasi-particle probability densities for the transition 4h 
5h where the tunneling hole has been forced to occupy the P2 
state. (Right) Result of the calculation where the hole tunnels into a 
_D-state. 



FIG. 5: Figure can be found at http://www.sst.nrel.gov/nano_pub/mts_preprint.pdfvColor online) Calculated quasi-particle probability 
densities for holes in dot D3 with a circular base. The reciprocal lattice vectors span a range of -7 to +1 10* m~^. 



